Estimation of the configuration of an optical element for an optical writing device

ABSTRACT

A method and an device for estimating the configuration of an optical element for an optical writing device are disclosed. The optical characteristics of the optical element are estimated on the basis of configuration data derived from actual measurement. The estimated optical characteristics are used to estimate the performance of the optical element. The optical element has a sufficiently small curve on the image surface of the writing device ascribable to the shift of focal point.

BACKGROUND OF THE INVENTION

[0001] The present invention relates to a method of estimating the configuration of an optical element for use in an optical writing device included in an image forming apparatus, an optical element, and a method of producing the same

[0002] Generally, a digital copier, laser printer or similar image forming apparatus includes an optical writing device that includes scanning optics. Today, it is a common practice with the scanning optics to use an aspherical lens represented by: $\begin{matrix} {{x(h)} = {\frac{{Ch}^{2}}{1 + \sqrt{1 - {\left( {1 + k} \right)C^{2}h^{2}}}} + {\sum{e_{i}h^{i}}}}} & {{Eq}.\quad (1)} \end{matrix}$

[0003] where h denotes a lens height, C denotes a paraxial curvature, k denotes a conic constant, and e_(i) denotes the coefficient of a polynomial.

[0004] An aspherical lens is mainly implemented as a plastic lens produced by injection molding.

[0005] To meet the increasing demand for high image quality and low cost scanning optics, it is necessary to reduce the number of optical elements constituting the scanning optics. This can be done with, e.g., a toric surface that is aspherical in the main scanning direction, but has a curvature varying in accordance with the main scanning direction. A lens with a toric surface is expressed as: $\begin{matrix} {{z\left( {x,y} \right)} = {\frac{{Cx}^{2}}{1 + \sqrt{1 - {\left( {1 + k} \right)C^{2}x^{2}}}} + {\sum\limits_{j = 0}^{m}{\sum\limits_{i = 0}^{n}{e_{ji}x^{i}y^{j}}}}}} & {{Eq}.\quad (2)} \end{matrix}$

[0006] where x denotes coordinates in the main scanning direction, y denotes coordinates in the subscanning direction, z denotes coordinates in the direction of an optical axis, c denotes a curvature in the main scanning direction, k denotes a conic constant, and e_(ji) denotes the coefficient of a polynomial.

[0007] An ultraprecision, free curved surface machining apparatus has made it possible to implement the opcical function surface of a lens as a toric surface with a high degree of freedom.

[0008] Each optical function surface of a lens or optical element included in the scanning optics usually has an effective range extending over several ten millimeters to several hundred millimeters in the main scanning direction, but only over several ten millimeters in the subscanning direction. When such a lens is produced by, e.g., plastic injection molding, a deviation from a designed value (configuration error hereinafter) occurs due to, e.g., the uneven contraction of resin.

[0009] Resins in general have a contraction ratio of about 0.7%. To fabricate a specular surface frame for injection molding in accordance with the designed configuration of an optical function surface, it has been customary to implement a configuration enlarged by similar enlargement using a contraction ratio that is determined either by theory or by experience. However, because an optical function surface is as great as several hundred millimeters, the nonlinear contraction of resin introduces a configuration error of several micrometers to several ten micrometers in the optical function surface, noticeably effecting the performance of the optics. In this respect, the estimation of the configuration error of a lens is essential. It is a common practice to measure the contour configuration of a molding in a plane that passes the design original in the mains scanning direction (generator) by use of a contour measuring device and then estimate a configuration error in accordance with the measured contour configuration. Typical of the contour measuring device is FORM TALYSURF available from Rank Taylor.

[0010] Japanese Patent Laid-Open Publication Nos. 6-129944 and 7-35541 each disclose a specifc estimation method for determining whether or not a configuration error is sufficiently small. One of conventional estimation methods optimizes the paraxial curvature radius of an aspherical equation in such a manner as to minimize the square sum of a configuration error (best fit R hereinafter). Another conventional method uses as an estimation parameter a difference between the best fit R and an aspherical equation, i.e., a configuration error. However, as the amount of asphyericality, i.e., a deviation from a designed spherical surface increases, the correlation between the belt fit R or the parameter for estimation and the optical performance decreases. Consequently, when a tolerance is distributed in order to prevent performance from being rejected by, e.g., optical simulation, some lenses may be rejected by configuration estimation, but may be allowed as to optical characteristic.

[0011] Japanese Patent Laid-Open Publication No. 9-89713, for example, proposes an estimation method using a parameter more closely correlated to optical characteristics than the above-discussed parameters. This estimation method produces the second-order derivative of a configuration error or a quadratic differential based on a difference between adjoining coordinates data or a difference between them. The method then confines the second-order derivative, the quadratic differential or the difference in a particular range so as to guarantee lens performance. The above document, however, simply teaches that the estimation method implements a closer correlation than the other conventional methods by presenting actual specific moldings. How the above method is effective for lens surfaces other than the actual moldings is not known.

[0012] On the other hand, Japanese Patent Laid-Open Publication Nos. 5-96572 and 7-60857 each propose to reduce the configuration error of a molding by measuring the configuration of the molding to thereby determine a configuration error and then correcting a specular surface frame in such a manner as to cancel the configuration error. This, however, requires a specular surface frame matching with all optical function surfaces to be corrected until the configuration error becomes sufficiently small.

[0013] Japanese Patent Laid-Open Publication Nos. 10-288749 and 11-77842 address to a refractive index distribution particular to a plastic lens produced by injection molding. Specifically, a non-uniform refractive index distribution in resin derives a lens effect and causes the focal point of optics to vary. To solve this problem, it is necessary to devise design or to execute an extra step for the reduction of the refractive index distribution after molding. Even if all optical elements are configured as designed, an optical scanning device including the optical elements cannot achieve expected optical performance unless the refractive index distribution is reduced to zero. This is because the shift of the focal point causes an image surface to curve. Particularly, assume that a spot diameter on an image surface is reduced to implement high density, scanning optics. Then, the F number of the optics must be reduced because the spot diameter is substantially proportional to the F number. This, however, causes the focal depth to decrease in proportion to the square of the spot diameter Therefore, the prerequisite with the development of such scanning optics is that the curve of the image surface ascribable to the shift of a focal point in the image surface be sufficiently small.

SUMMARY OF THE INVENTION

[0014] It is therefore a first object of the present invention to provide a configuration estimating method capable of estimating the optical characteristics of a lens on the basis of the result of measurement and estimating the performance of the lens on the basis of the estimated characteristics and therefore adaptive to any desired lens surface and closely correlated to the optical characteristics, and a device for practicing the same.

[0015] It is a second object of the present invention to provide a method of producing an optical element capable of sufficiently reducing the curve of an image surface ascribable to the shift of a focal point on the image surface and producing such an optical element in a short period of time at low cost, and the optical element.

[0016] In accordance with the present invention, a method of estimating the surface configuration of a lens begins with a step of measuring the contour configuration of the lens. The configuration error of the lens that is a deviation from a designed configuration is determined. Subsequently, a partial curvature at each lens height is determined on the basis of the configuration error. Thereafter, a curvature proportional coefficient, which is a shift of a focal point for a unit curvature at each lens height, as measured in the direction of an optical axis on an image surface, is determined. Finally, the shift is estimated on the basis of the curvature proportional coefficient and partial curvature.

BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The above and other objects, features and advantages of the present invention will become more apparent from the following detailed description taken with the accompanying drawings in which:

[0018]FIG. 1 is a view showing an, optical writing unit included in an image forming apparatus and to which the present invention is applicable;

[0019]FIG. 2 is a view showing optics adjoining a luminous flux corresponding to the rotation angle of a polygon motor;

[0020]FIG. 3 is a flowchart demonstrating a specific procedure for estimating the shift of a focal point unique to a first embodiment of the present invention;

[0021]FIG. 4 is a flowchart demonstrating another specific procedure for estimating the shift of a focal point available with the first embodiment of the present invention; and

[0022]FIG. 5 is a flowchart showing a specific procedure for reducing the curve of an image surface unique to a second embodiment of the present invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0023] Referring to FIG. 1 of the drawings, an optical writing device to which the present invention is applied is shown. As shown, the optical writing device includes a semiconductor laser or similar light source 1. A luminous flux issuing from the light source 1 is transmitted through a collimator lens 2 to become a parallel flux. The parallel flux is then converged by a cylindrical lens 3 onto the reflection face of a polygonal mirror 4 in the form of a line. The polygonal mirror 4 in rotation deflects the incident flux toward an image surface 6 via scanning optics 5. A photoconductive element, not shown, has its surface positioned at the image surface 6.

[0024] Assume that the direction in which the polygonal mirror 4 deflects the luminous flux is a main scanning direction, and that the direction perpendicular to the main scanning direction is a subscanning direction. Further, assume that a position on each optical function surface in the main scanning direction is a lens height, and that a position on the image surface 6 in the main scanning direction is an image height.

[0025] Preferred embodiments of the present invention will be described hereinafter with reference to the accompanying drawings.

First Embodiment

[0026] Generally, when the configuration of a lens is close to a designed configuration, the shift of a focal point on an image surface in the direction of optical axis is substantially proportional to the partial curvature of the configuration error of a lens surface at each lens height. By using this relation, the illustrative embodiment measures the contour configuration of a subject lens, determines a deviation of the measured configuration from a designed configuration, then determines a partial curvature of the deviation, then estimates the shift of a focal point at any des red lens height, and then evaluates the lens.

[0027] First, assume that all lenses have ideal configurations. Then, at a given time when a luminous flux passes through part of the effective range of a given lens at a particular position of the lens in accordance with the rotation angle of a polygon motor, the flux is focused on a fixed position adjoining an image surface. This may be interpreted to mean that optical elements are independent of each other at a given time. More specifically, it is possible to determine, e.g., the focal distance of the individual optical element. In addition, because the individual lens has a smoothly contiguous surface, the focal distance and other optical characteristics are considered to be also smoothly contiguous. Although each lens may have an ideal configuration, the focal point at a given time is, in a strict sense, shifted from the image surface, as will be discussed specifically later. At this stage of description, consideration is given to a shift from a focal point to occur in the ideal condition when a subtle configuration error is introduced in a lens having an ideal configuration.

[0028] Second, assume that a lens to be estimated involves a configuration error. A configuration error is sufficiently small at any time with respect to a designed configuration. Therefore, assume that the influence of the configuration error in a luminous flux on optical characteristics at any time lies in a range that can be approximated by the paraxial theory. Today, the shift of a focal point ascribable to a configuration error is a decisive factor that effects yield. The shift of a focal point is ascribable to variation in partial curvature, as confirmed by optical simulations. In light of this, the illustrative embodiment estimates the shift of a focal point on the image surface ascribable to variation in partial curvature.

[0029] The influence of a configuration error on the front and rear of a given lens will be described by using Snell laws. A Snell's equation in an ideal condition is: $\begin{matrix} {\frac{n^{\prime}}{s_{d}^{\prime}} = {\frac{n}{s} + {\left( {n^{\prime} - n} \right)C_{d}}}} & {{Eq}.\quad (3)} \end{matrix}$

[0030] where s denotes an imaging position at the object side, s_(d)′ denotes an imaging position at the image side, n denotes a refractive index at the object side, and C_(d) denotes the curvature of a lens surface.

[0031] As shown in FIG. 2, assume that the curvature C_(d) is varied to C_(r) by a configuration error. Then, the resulting focal point s_(t)′ at the image side is expressed as: $\begin{matrix} {\frac{n^{\prime}}{s_{r}^{\prime}} = {\frac{n}{s} + {\left( {n^{\prime} - n} \right)C_{r}}}} & {{Eq}.\quad (4)} \end{matrix}$

[0032] The shift Δ_(s)′ of the focal point at the image side is produced by: $\begin{matrix} \begin{matrix} {{\Delta \quad s^{\prime}} = \quad {{s_{r}^{\prime} - s_{d}^{\prime}} = {{- s_{r}^{\prime}}s_{d}^{\prime}\frac{n^{\prime} - n}{n^{\prime}}\left( {C_{r} - C_{d}} \right)}}} \\ {\approx \quad {{- s_{d}^{\prime 2}}\frac{n^{\prime} - n}{n^{\prime}}\Delta \quad C}} \end{matrix} & {{Eq}.\quad (5)} \end{matrix}$

[0033] Because a curvature error ΔC is sufficiently smaller than the designed curvature C_(d), the shift of the focal point Δ_(s)′ is smaller than the focal distance s′ and approximately expressed as:

S_(r) ¹≈S_(d) ¹   Eq. (6)

[0034] So long as lenses following the above lens surface have designed values, a magnification is constant, so that the shift of the focal point on the image surface is proportional to the curvature error. In this case, a proportional coefficient is a constant dependent on the designed configuration.

[0035] It is therefore possible to estimate the shift of a focal point of a given lens on the image surface by producing a proportional coefficient, which is the shift of a focal point on the image surface for a unit curvature at any lens height, and a curvature error at any lens height. FIG. 3 demonstrates a procedure for dealing with measured data.

[0036] The result of measurement is output in the form of dot-sequential data. That is, each dot is output as coordinates data in a coordinates system that is dependent on a measuring apparatus. To, separate a designed configuration component from the measured data, the coordinates system used for measurement and the coordinates system used for design must be identical. In practice, however, such two systems cannot be identical due to, e.g., the accuracy of mounting of a subject lens. Coordinates conversion is therefore necessary for representing the dot sequence of the coordinate system used for measurement by the coordinates system used for design. Because a matrix for coordinates conversion cannot be directly produced, the coordinates system used for measurement is replaced with a coordinates system whose error is smallest. This can be done by optimizing the above matrix in such a manner as to minimize the configuration error. As for bidimensional data, among three degrees of freedom, i.e., a shift sx in the direction of optical axis, a shift sy in the direction perpendicular to the direction of optical axis and a rotation angle θ, at least sy and θ are optimized. In practice, to correct the mounting error, tridimensional optimization may be effected by using the square sum of the configuration error after coordinates conversion as an estimation function.

[0037] A partial curvature is produced from the configuration error, as follows. First, the configuration error data is divided into a plurality of dot sequences. Each divided dot sequence has a length in the y axis direction that should preferably be close to the width of the luminous flux. Subsequently, the divided dot sequences each are approximated to a spherical surface by a minimum square scheme in order to determine the curvature of the spherical surface. For the approximation, use is made of an equation representative of a spherical surface: $\begin{matrix} {{x - {x0}} = \frac{{C\left( {y - {y0}} \right)}^{2}}{1 + \sqrt{1 - {C^{2}\left( {y - {y0}} \right)}^{2}}}} & {{Eq}.\quad (7)} \end{matrix}$

[0038] where (x0, y0) denotes the coordinates of the apex of the spherical surface, and c denotes the curvature of the spherical surface. In practice, the curvature c of the spherical surface is produced by optimizing such parameters by use of the minimum square scheme.

[0039] Another specific method of producing a partial curvature from the configuration error is as follows. First, the configuration error data are divided into a plurality of dot sequences. Again, each divided dot sequence has a length in the y axis direction that should preferably be close to the width of the luminous flux. Subsequently, the divided dot sequences each are approximated to a quadratic polynomial by the minimum square scheme. The square term of the quadratic polynomial is doubled to produce a curvature. In practice, use is made of a singular value resolving method or similar linear minimum square method.

[0040] Still another specific method of producing a partial curvature from the configuration error is as follows. First, the configuration error data are approximated to a polynomial (8), which is shown below, having a suitable degree so as to produce the coefficients of the individual terms. Subsequently, the coefficients of a second-order derivative (9), which is also shown below, of the approximated polynomial are produced. $\begin{matrix} {{f(y)} = {\sum\limits_{i = 0}^{n}{e_{i}y^{i}}}} & {{Eq}.\quad (8)} \\ {{f^{''}(y)} = {\sum\limits_{i = 2}^{n}{{i\left( {i - 1} \right)}e_{i}y^{i - 2}}}} & {{Eq}.\quad (9)} \end{matrix}$

[0041] The above operation is repeated with a plurality of lens heights.

[0042] A specific method of estimating the shift of a focal point will he described hereinafter. For preparation, the sample configuration error of a suitable model equation is added to the designed value of a particular lens surface. Optical simulation is then effected by using a plurality of deflector angles. Subsequently, a lens height of a point where a chief ray derived from each deflector angle and the lens surface meet each other, as measured in the main scanning direction, and a shift from the designed focal point on, the image surface are produced. Also, the partial curvature of a polynomial configuration error at each lens height is determined by the previously stated method. The shift of the focal point at each lens height is divided by the partial curvature to thereby produce a partial curvature proportional coefficient, which is the shift of the focal point for a unit curvature. Thereafter, the relation between the lens height and the partial curvature proportional coefficient is approximated to a suitable function, e.g., a polynomial. By storing the parameter of the polynomial, it is possible to determine the proportional coefficient of a shift at any lens height.

[0043] To estimate the shift of a focal point, a partial curvature and a proportional coefficient are produced at each of a plurality of lens heights from the results of measurement of actual samples and then added. FIG. 4 demonstrates the procedure for dealing with measured data.

[0044] Another specific method of estimating the shift of a focal point is as follows. Again, for preparation, the sample configuration error of a suitable model equation is added to the designed value of a particular lens surface. Optical simulation is then effected by using a plurality of deflector angles. Subsequently, there are produced a lens height at a point where a chief ray derived from each deflector angle and the lens surface meet each other, as measured in the main scanning direction, an image height on the image surface and a shift from the designed focal point on the image surface. Also, the partial curvature of a polynomial configuration error at each lens height is determined by the previously stated method. Subsequently, the relation between the lens height and the image height is approximated to a model equation, e.g., a polynomial, and a parameter thereof is preserved. Further, the shift of focal point at each lens height is divided by the partial curvature in order to produce a proportional coefficient, which is the shift of a focal point for a unit curvature. Such proportional coefficients are approximated to a suitable function, e.g., a polynomial, and a parameter thereof is preserved. This allows the proportional coefficient of the shift of a focal point at any lens height to be produced. Thereafter, the product of the image height and partial curvature and the proportional coefficient is produced at each of a plurality of lens heights from the result of actual measurement, thereby determining a relation between the image height and the shift of focal point.

[0045] When the sample configuration error is spherical and has a certain curvature, the partial curvature of the configuration error at any lens height is identical. In this case, a partial curvature does not have to be produced. That is, the shift of the focal point is divided by a curvature or constant to thereby produce a partial curvature.

[0046] When the sample configuration error is parabolic, the partial curvature is constant at all lens heights if use is made of the previously stated method that produces a second-order derivative. It is therefore not necessary to determine the partial curvature of the sample configuration error as when the sample configuration error is spherical.

[0047] Assume that the shift of the focal point is estimated in terms of the relation between the image height in the image surface and the shift. Then, the influence of the first surface and that of the second surface of a lens are estimated as image surface curves, which are the shifts of the focal point in a single image surface. By adding the image surface curvatures, it is possible to estimate the shift of the focal point of a single lens.

[0048] If a personal computer, for example, capable of executing the above procedure is connected to the measuring apparatus, then the procedure can be implemented by software on the personal computer.

[0049] As stated above, the illustrative embodiment has various unprecedented advantages, as enumerated below.

[0050] (1) A shift of the focal point can be estimated from a configuration error, so that a lens can be evaluated on the basis of the shift. Consequently, the result of configuration evaluation and the optical performance of a lens are more closely correlated to each other. This increases yield, allows the tolerances of various parameters for machining to be increased, and reduces the production cost of lenses.

[0051] (2) Errors unavoidable when a subject lens is mounted to a jig can be corrected to promote evaluation with high repeatability.

[0052] (3) Measurement data derived from a contour configuration are divided into partial data each having a length substantially equal to the diameter of a luminous flux. A curvature is then determined with each of the partial data. It is therefore possible to more accurately estimate a configuration as to a local configuration error existing in cart of a lens surface.

[0053] (4) Approximation to a quadratic polynomial is more stable and implements estimation in a shorter time than approximation to a spherical surface.

[0054] (5) A polynomial promotes more stable configuration estimation by reducing the influence of, e.g., fine dust present on the surface of a subject lens during measurement. Dust and other deposits present on a lens would cause coordinates to locally vary to a noticeable degree.

[0055] (6) The proportional coefficient of a shift of the focal point on an image surface for a unit curvature is determined by optical simulations beforehand. A relation between the proportional coefficient and the lens height can be approximated to a suitable function and stored. This allows a proportional coefficient to be easily produced for any desired lens height and therefore realizes rapid estimation.

[0056] (7) A relation between the image height on the image surface and a shift of the focal point is finally obtained. By comparing such a relation with the result of measurement effected with an actual lens, it is possible to examine the accuracy of estimation and to correct the proportional coefficient.

[0057] (8) Because sample configuration errors are implemented by a spherical surface, a procedure for estimating the curvatures of a sample configuration is simplified.

[0058] (9) Because sample configuration errors are implemented by a parabolic surface, a procedure for estimating the curvatures of a sample configuration is simplified. Further, while the quadratic term of a polynomial is included in the non-spherical coefficients of an optical simulator, it is, in many cases, not used and allows a numerical value to be input. This makes it needless to modify a simulator or needs a minimum of modification.

[0059] (10) By adding a shift of the focal point of the first surface and that of the second surface of a lens, it is possible to estimate the total accuracy of the lens. For example, even a lens whose individual surface does not satisfy a tolerance, but satisfies it as to total accuracy, can be adequately estimated. This increases yield.

[0060] (11) A configuration estimating device with a close correlation to optical characteristics is achievable.

Second Embodiment

[0061] This embodiment executes the procedure of the previous embodiment in the reverse way. Briefly, the illustrative embodiment measures shifts of the focal point at different image heights, generates, based on the measured shifts, a corrected configuration for canceling the shift in a particular optical function surface, and correcting a specular surface frame for molding to thereby reduce the field curvature of the entire optics. The illustrative embodiment will be described specifically with reference to FIG. 5.

[0062] Initial Molding

[0063] First, a mold is produced by a conventional method and used to produce optical elements by injection molding. At this instant, molding conditions are accurately examined to stably produce optical elements in the same configuration.

[0064] Measurement of Shift of Focal Point

[0065] The optical elements are built in an optical scanning device in order to measure shifts of the focal point. Specifically, the polygonal mirror is fixed in place while a device for measuring a beam to issue from a laser is mounted on, e.g., a direct-acting stage. The stage is moved in the direction of optical axis in order to measure the spot diameter of a luminous flux. The above device measures the beam profile of a luminous flux, incident to a particular surface and outputs, e.g., spot diameters in the main and subscanning directions. This produces deviations from the image surface at various image heights and the distributions of spot diameters in the main and subscanning directions. These factors are approximated to, e.g., a quadratic function to thereby produce minimum values, each of which is representative of a shift of the focal point at a particular image height. Consequently, there are produced an image height distribution and shift distributions in the main and subscanning directions.

[0066] Calculation of Lens Height

[0067] Lens heights corresponding to the image heights, which are used to determine the shifts of the focal point, are calculated on a particular optical function surface. Specifically, a distribution of image heights and that of lens heights are determined by simulations beforehand. Such distributions are approximated to a function, and the coefficient of the function is stored. The function is then solved to determine lens heights.

[0068] Calculation of Sensitivity

[0069] There are calculated, on a particular optical function surface, the sensitivity of partial curvatures corresponding to the image heights from which the shifts of the focal point are produced. Specifically, a distribution of image heights and that of sensitivity of partial curvatures are determined by simulation beforehand. More specifically, the following procedure is executed by varying the deflection angle while giving a configuration error only to a particular surface. An image height at which the chief ray of the luminous flux is incident to the image surface is determined. Configuration errors in a range which the chief ray passes through a particular optical function surface are approximated to a spherical surface or a quadratic curve, thereby determined a partial curvature selected. Spot diameters in the main and subscanning directions are produced at a plurality of positions close to the image surface. Such factors each are approximated to, e.g., a particular quadratic function to thereby determine a minimum value, which is representative of a shift of the focal point. The shifts in the main and subscanning directions each are divided by associated one of the partial curvatures for thereby determining sensitivity, which is a shift of the focal point for a unit partial curvature. The sensitivity is approximated to a function, and the coefficient of the function is stored. In practice, the individual function is solved.

[0070] Calculation of Partial Curvature

[0071] The sensitivity and the shift of the focal point are used to produce a corrected partial curvature that cancels the shift of the focal point of the entire optics in a particular optical function surface. A corrected partial curvature dC is produced by: $\begin{matrix} {{dC} = {- \frac{d}{k}}} & {{Eq}.\quad (10)} \end{matrix}$

[0072] where d and k denote a shift of the focal point and sensitivity, respectively.

[0073] Calculation of Corrected Configuration

[0074] Finally, lens heights in one or both of the main and subscanning directions and the distribution of corrected partial curvatures are used to calculate a corrected configuration in the particular optical function surface. For this purpose, the coefficient of a model equation representative of a curved surface is determined, or alternatively map data in the optical function surface are generated. A corrected configuration is output such that a height at any position can be calculated for producing NC (Numerical Control) data.

[0075] A specific method of producing a corrected configuration is as follows. A surface on which a partial curvature in the subscanning direction varies in the main scanning direction is produced by: $\begin{matrix} {C_{s} = {f(x)}} & {{Eq}.\quad (11)} \\ {{z\left( {x,y} \right)} = \frac{C_{s}y^{2}}{1 + \sqrt{1 - {C_{s}^{2}y^{2}}}}} & {{Eq}.\quad (12)} \end{matrix}$

[0076] It follows that if the partial curvature distribution in the subscanning direction can be approximated to a function f(x), then a corrected configuration in the subscanning direction can be directly produced. As for the function, use should preferably be made of a polynomial or a spline.

[0077] Another specific method of producing a corrected configuration is as follows. Generally, when a partial curvature is small, it is substantially equal to the two-step partial differential of a configuration. By using this fact, the specific method to be described integrates the partial curvature distribution for thereby producing a corrected configuration. Specifically assume that the function model of a corrected configuration to be produced is the following n×2 degree polynomial: $\begin{matrix} {{z\left( {x,y} \right)} = {\sum\limits_{j = 0}^{2}{\sum\limits_{i = 0}^{n}{e_{ji}x^{i}y^{j}}}}} & {{Eq}.\quad (13)} \end{matrix}$

[0078] Two-step partial differential is expressed as: $\begin{matrix} {\frac{\partial^{2}z}{\partial x^{2}} = {\sum\limits_{j = 0}^{2}{\sum\limits_{i = 2}^{n}{{i\left( {i - 1} \right)}e_{ji}x^{i - 2}y^{j}}}}} & {{Eq}.\quad (14)} \\ {\frac{\partial^{2}z}{\partial y^{2}} = {{\sum\limits_{j = 2}^{2}{\sum\limits_{i = 2}^{n}{{j\left( {j - 1} \right)}e_{ji}x^{i}y^{j - 2}}}} = {\sum\limits_{i = 0}^{n}{2e_{2i}x^{i}}}}} & {{Eq}.\quad (15)} \end{matrix}$

[0079] It is therefore possible to produce a corrected configuration by producing e_(ij) from partial curvatures in the main and subscanning directions. As for the subscanning direction, the lens height distribution and partial curvature distribution are approximated to n-degree polynomials, and the coefficients of the polynomials each are multiplied by 0.5 to thereby produce respective e_(2i). As for the main scanning direction, the lens height distribution and partial curvature distribution each are approximated to an n-2-degree polynomial and then integrated two times to thereby produce respective e_(0i). Assume that NC data for machining a specular surface frame for correction should be generated. Then, the above procedure allows a height at any desired height to be produced and therefore allows a corrected configuration in a particular optical function surface that cancels the configuration errors of the entire optics to be produced.

[0080] A further specific method of producing a corrected configuration is as follows. This method calculates a corrected configuration on the basis of the partial curvature distribution in the main scanning direction. First, the lens height distribution and partial curvature distribution each are numerically integrated two times, thereby producing a lens height distribution and a corrected configuration distribution. Such distributions each are approximated to a suitable function model. This allows any desired function model other than a spline, trigonometric function or similar polynomial model included in design equations to be used.

[0081] Because a corrected configuration is produced by integration, the influence of, e.g., focal point measurement errors often accumulate. In light of this, when the entire corrected configuration has an excessively great size, it is preferable to correct low-degree coefficients, which influence partial curvature little, for thereby reducing the absolute value of the corrected configuration.

[0082] Machining for Correction

[0083] The corrected configuration produced by the procedure described above is superposed on a designed configuration so as to produce a new configuration of the optical function surface. The specular surface frame of a mold is produced on the basis of the new configuration. Assume that use is made of a free curved surface machining apparatus for producing an absolute configuration by, e.g., cutting or grinding. Then, after the corrected configuration has been superposed on the designed configuration, NC data for machining a specular surface frame by a conventional method and used for machining. On the other hand, assume that use is made of a grinding machine capable of removing a work by a specified relative amount, the corrected configuration is directly subjected to similar transformation and plus-minus inversion in the direction of height to thereby generate NC data for machining.

[0084] Molding

[0085] Finally, injection molding is effected in the same conditions as the initial molding. As a result, despite that all the optical function surfaces of an optical element are deviated from designed configurations, the shift of focal point of the optical element, except for a scatter component, is substantially zero when the element is mounted to an optical scanning device.

[0086] The illustrative embodiment will be described more specifically hereinafter.

[0087] (1) The above-described procedure cannot reduce the shift of the focal point to zero when executed only once because of the machining errors of a specular surface frame and errors ascribable to molding. In light of this, the measurement of a shift stated earlier is executed with the produced optical element to see if it satisfies target performance required of an optical scanning device. If the performance of the optical element is short of the target performance, then the procedure is repeated until the former reached the latter.

[0088] (2) A particular optical function surface should preferably be a surface whose sensitivity is generally low in both of the main and subscanning directions. This is because the corrected configuration of such a surface is large sized and therefore little susceptible to the influence of the absolute accuracy of a machining apparatus. In many cases, scanning optics does not include an optical function surface whose sensitivity is low in both of main and subscanning directions. In such a case, a particular surface in each of the main and subscanning directions, i.e., two surfaces should preferably be corrected.

[0089] Further, when it is desired to correct only one specular surface frame, at least one surface in the scanning optics may be made flat. Particularly, it is preferable to additionally design a parallel flat sheet and mold an optical element therewith by injection molding. To prevent the sensitivity of the resulting molding from increasing in sensitivity in the main and subscanning directions, the molding should preferably be positioned close to the image surface. In this condition, a corrected configuration is produced in accordance with the above flat surface, and then the specular surface frame thereof is machined for correction. Machining a flat surface for correction by several micrometers to several ten micrometers is easier and more accurate than machining an optical function surface, which is curved. Moreover, use can be made of a machining apparatus with a small stroke that is far lower in cost than an ultraprecision, free curved surface machining apparatus. In this manner, only if the substantially flat specular surface is corrected for a single optical function surface, the shift of focal point of a scanning device can be more accurately corrected.

[0090] (3) The procedure beginning with the calculation of lens heights and ending with the calculation of a corrected configuration can be implemented by, e.g., software installed in a personal computer. It is therefore possible to produce software that sequentially executes the procedure shown and described and outputting NC data for correction by inputting measured shifts of focal point.

[0091] As described above, the illustrative embodiment achieves various unprecedented advantages, as enumerated below.

[0092] (1) When plastic optical elements unique to the illustrative embodiment are built in an optical scanning device, the shift of the focal point except for scatter is reduced to substantially zero. The scanning device therefore achieves high performance. In addition, machining for correction is effected on the basis of the shifts of the focal point of the scanning device actually measured, the influence of a refractive index distribution is corrected at the same time. It follows that a refractive index distribution, which may be left in the optical elements, does not effect the performance of the scanning device.

[0093] (2) Shifts of the focal point are actually measured and then corrected by use of optical function surfaces that have margins as to machining accuracy. Therefore, even if the machining accuracy available is short of accuracy required of a particular optical function surface of the optics, optical elements particularly feasible for a scanning device whose beam diameter is small can be produced. This is desirable from the industrial application standpoint.

[0094] (3) The total accuracy of a lens is estimated as well before the correction of a specular surface frame. Therefore, even if a single surface of a lens does not satisfy a tolerance, it can be so corrected to satisfy the tolerance as to total accuracy and therefore achieves high accuracy.

[0095] (4) Optical elements for an optical scanning device that achieve target performance can be produced.

[0096] (5) The procedure ends in a short time because optical simulation does not have to be conducted each time.

[0097] (6) The configuration of an optical function surface that cancels the shift of the focal point of optics is attainable without resorting to additional optical simulations or the additional design of optical function surfaces.

[0098] (7) NC data is attainable in a short period because only part of the coefficients of design equations should be changed for the calculation of a corrected configuration and the generation of NC data.

[0099] (8) The curve of an image surface can be corrected without resorting to an ultraprecision, free curved surface machining device having a sufficient operation range. This successfully reduces the machining cost of a specular surface frame.

[0100] (9) NC data for correction machining can be almost automatically produced from data representative of measured shifts of the focal point. This reduces designer's errors and allows correction to be completed in a short period.

[0101] Various modifications will become possible for those skilled in the art after receiving the teachings of the present disclosure without departing from the scope thereof. 

What is claimed is:
 1. A method of estimating a surface configuration of a lens, comprising: a configuration measuring step of measuring a contour configuration of the lens; an error determining step of determining a configuration error of the lens that is a deviation from a designed configuration; a partial curvature determining step of determining a partial curvature at each lens height on the basis of the configuration error; and a shift estimating step for determining a curvature proportional coefficient, which is a shift of a focal point for a unit curvature at each lens height, as measured in a direction of an optical axis on an image surface, and estimating said shift on the basis of said curvature proportional coefficient and the partial curvature.
 2. The method as claimed in claim 1, wherein said error determining step comprises: an error correcting step of determining, after translational and rotational coordinates conversion with all or part of coordinates data, a deviation from the designed configuration, and determining optimal values of part of or all of translation amounts and rotation amounts that minimize a square sum of said deviation; and a designed configuration separating step of determining a deviation of a configuration subjected to optimal translation and optimal rotation from the designed configuration.
 3. The method as claimed in claim 2, wherein said partial curvature determining step comprises dividing configuration error data, and determining a curvature by approximating each divided configuration error data to a spherical surface by a minimum square method.
 4. The method as claimed in claim 2, wherein said partial curvature determining step comprises dividing configuration error data, approximating each divided configuration error data to a polynomial by a minimum square method, and doubling a square term of said polynomial to thereby produce a partial curvature.
 5. The method as claimed in claim 2, wherein said partial curvature determining step comprises approximating configuration error data to a polynomial, and producing a second-order derivative of said polynomial.
 6. The method as claimed in claim 2, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of a focal point in a direction of an optical axis in a corresponding image surface, determining partial curvatures of said sample configuration errors, determining partial curvature proportional coefficients each being representative of a shift of the focal point for a unit curvature and corresponding to a particular lens height, and approximating a resulting relation to a function while storing a parameter of said function; and a shift producing step of solving an approximate function by using the parameter to thereby produce partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 7. The method as claimed in claim 6, wherein the sample configuration errors are spherical.
 8. The method as claimed in claim 6, wherein the sample configuration errors are parabolic.
 9. The method as claimed in claim 2, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights and corresponding image heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of the focal point in a direction of an optical axis in a corresponding image surface, approximating said lens heights and said image heights to a function while storing a parameter of said function, determining partial curvatures of said sampled configuration errors, determining partial curvature proportional coefficients each being representative of a shift of a focal point for a unit curvature and corresponding to a particular lens height, and approximating a relation between said partial curvature proportional coefficients and said lens heights while storing a parameter of said function; and a shift producing step of solving approximate functions by using the parameter to thereby produce image heights and partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 10. The method as claimed in claim 1, wherein said partial curvature determining step comprises dividing configuration error data, and determining a curvature by approximating each divided configuration error data to a spherical surface by a minimum square method.
 11. The method as claimed in claim 10, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of a focal point in a direction of an optical axis in a corresponding image surface, determining partial curvatures of said sample configuration errors, determining partial curvature proportional coefficients each being representative of a shift of the focal point for a unit curvature and corresponding to a particular lens height, and approximating a resulting relation to a function while storing a parameter of said function; and a shift producing step of solving an approximate function by using the parameter to thereby produce partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 12. The method as claimed in claim 11, wherein the sample configuration errors are spherical.
 13. The method as claimed in claim 11, wherein the sample configuration errors are parabolic.
 14. The method as claimed in claim 10, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights and corresponding image heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of the focal point in a direction of an optical axis in a corresponding image surface, approximating said lens heights and said image heights to a function while storing a parameter of said function, determining partial curvatures of said sampled configuration errors, determining partial curvature proportional coefficients each being representative of a shift of a focal point for a unit curvature and corresponding to a particular lens height, and approximating a relation between said partial curvature proportional coefficients and said lens heights while storing a parameter of said function; and a shift producing step of solving approximate functions by using the parameter to thereby produce image heights and partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 15. The method as claimed in claim 10, further comprising: a step of approximating a relation between an image height of a first surface of a subject lens and a shift of the focal point to a model equation; a step of approximating a relation between an image height of a second surface of the subject lens by a model equation; and a step of producing a sum of said model equations to thereby determine a total shift of the focal point of the subject lens.
 16. The method as claimed in claim 1, wherein said partial curvature determining step comprises dividing configuration error data, approximating each divided configuration error data to a polynomial by a minimum square method, and doubling a square term of said polynomial to thereby produce a partial curvature.
 17. The method as claimed in claim 16, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of a focal point in a direction of an optical axis in a corresponding image surface, determining partial curvatures of said sample configuration errors, determining partial curvature proportional coefficients each being representative of a shift of the focal point for a unit curvature and corresponding to a particular lens height, and approximating a resulting relation to a function while storing a parameter of said function; and a shift producing step of solving an approximate function by using the parameter to thereby produce partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 18. The method as claimed in claim 16, wherein the sample configuration errors are spherical.
 19. The method as claimed in claim 16, wherein the sample configuration errors are parabolic.
 20. The method as claimed in claim 16, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights and corresponding image heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of the focal point in a direction of an optical axis in a corresponding image surface, approximating said lens heights and said image heights to a function while storing a parameter of said function, determining partial curvatures of said sampled configuration errors, determining partial curvature proportional coefficients each being representative of a shift of a focal point for a unit curvature and corresponding to a particular lens height, and approximating a relation between said partial curvature proportional coefficients and said lens heights while storing a parameter of said function; and a shift producing step of solving approximate functions by using the parameter to thereby produce image heights and partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 21. The method as claimed in claim 16, further comprising: a step of approximating a relation between an image height of a first surface of a subject lens and a shift of the focal point to a model equation; a step of approximating a relation between an image height of a second surface of the subject lens by a model equation; and a step of producing a sum of said model equations to thereby determine a total shift of the focal point of the subject lens.
 22. The method as claimed in claim 1, wherein said partial curvature determining step comprises approximating configuration error data to a polynomial, and producing a second-order derivative of said polynomial.
 23. The method as claimed in claim 22, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of a focal point in a direction of an optical axis in a corresponding image surface, determining partial curvatures of said sample configuration errors, determining partial curvature proportional coefficients each being representative of a shirt of the focal point for a unit curvature and corresponding to a particular lens height, and approximating a resulting relation to a function while storing a parameter of said function; and a shift producing step of solving an approximate function by using the parameter to thereby produce partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 24. The method as claimed in claim 22, wherein the sample configuration errors are spherical.
 25. The method as claimed in claim 23, wherein the sample configuration errors are parabolic.
 26. The method as claimed in claim 22, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights and corresponding image heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of the focal point in a direction of an optical axis in a corresponding image surface, approximating said lens heights and said image heights to a function while storing a parameter of said function, determining partial curvatures of said sampled configuration errors, determining partial curvature proportional coefficients each being representative of a shift of a focal point for a unit curvature and corresponding to a particular lens height, and approximating a relation between said partial curvature proportional coefficients and said lens heights while storing a parameter of said function; and a shift producing step of solving approximate functions by using the parameter to thereby produce image heights and partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 27. The method as claimed in claim 22, further comprising: a step of approximating a relation between an image height of a first surface of a subject lens and a shift of the focal point to a model equation; a step of approximating a relation between an image height of a second surface of the subject lens by a model equation; and a step of producing a sum of said model equations to thereby determine a total shift of the focal point of the subject lens.
 28. The method as claimed in claim 1, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of a focal point in a direction or an optical axis in a corresponding image surface, determining partial curvatures of said sample configuration errors, determining partial curvature proportional coefficients each being representative of a shift of the focal point for a unit curvature and corresponding to a particular lens height, and approximating a resulting relation to a function while storing a parameter of said function; and a shift producing step of solving an approximate function by using the parameter to thereby produce partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 29. The method as claimed in claim 28, wherein the sample configuration errors are spherical.
 30. The method as claimed in claim 28, wherein the sample configuration errors are parabolic.
 31. The method as claimed in claim 28, further comprising: a step of approximating a relation between an image height of a first surface of a subject lens and a shift of the focal point to a model equation; a step of approximating a relation between an image height of a second surface of the subject lens by a model equation; and a step of producing a sum of said model equations to thereby determine a total shift of the focal point of the subject lens.
 32. The method as claimed in claim 1, wherein said shift estimating step comprises: a preparatory step of determining by an optical simulation lens heights and corresponding image heights when sample configuration errors of only a subject lens surface are superposed on designed values of a lens surface, and shifts of the focal point in a direction of an optical axis in a corresponding image surface, approximating said lens heights and said image heights to a function while storing a parameter of said function, determining partial curvatures of said sampled configuration errors, determining partial curvature proportional coefficients each being representative of a shift of a focal point for a unit curvature and corresponding to a particular lens height, and approximating a relation between said partial curvature proportional coefficients and said lens heights while storing a parameter of said function; and a shift producing step of solving approximate functions by using the parameter to thereby produce image heights and partial curvature proportional coefficients, and producing a shift curve on the image surface by using said partial curvature proportional coefficients.
 33. The method as claimed in claim 32, wherein the sample configuration errors are spherical.
 34. The method as claimed in claim 32, wherein the sample configuration errors are parabolic.
 35. The method as claimed in claim 32, further comprising: a step of approximating a relation between an image height of a first surface of a subject lens and a shift of the focal point to a model equation; a step of approximating a relation between an image height of a second surface of the subject lens by a model equation; and a step of producing a sum of said model equations to thereby determine a total shift of the focal point of the subject lens.
 36. The method as claimed in claim 1, further comprising: a step of approximating a relation between an image height of a first surface of a subject lens and a shift of the focal point to a model equation; a step of approximating a relation between an image height of a second surface of the subject lens by a model equation; and a step of producing a sum of said model equations to thereby determine a total shift of the focal point of the subject lens.
 37. A method of producing an optical element by molding, comprising: a shift measuring step of mounting the optical element molded by injection molding to an optical scanning device, which is equivalent to an optical scanning device to be actually used, and measuring shifts of a focal point in a direction of optical axis at a plurality of image heights on an image surface; a lens height calculating step of calculating positions on a particular optical function surface that correspond one-to-one to the plurality of image heights; a sensitivity calculating step of calculating sensitivity in terms of proportional coefficients of partial curvatures and shifts of the focal point on the particular optical function surface that correspond one-to-one to the plurality of image heights; a partial curvature calculating step of calculating partial curvatures on the particular optical function by using the shifts of the focal point and the sensitivity; a corrected configuration calculating step of outputting, based on the lens heights and the partial curvatures, a corrected configuration of the particular function surface in a form of coefficients of a function model or map data; a correction machining step of machining for correction a specular surface frame of a mold corresponding to the particular optical function surface; and a molding step of molding the optical element by using the specular surface frame corrected.
 38. The method as claimed in claim 37, further comprising a step of executing said shift measuring step by using the optical element corrected and repeating, if said optical element does not have a target configuration, said lens height calculating step through said molding step.
 39. The method as claimed in claim 38, wherein said lens height calculating step comprises determining by an optical simulation lens heights in a main scanning direction, at each of which a chief ray deflected by a particular angle is incident to a particular optical function surface, and image heights in said main scanning direction at each of which said chief ray is incident to an image surface, approximating a lens height distribution and an image height distribution to a polynomial or similar model while storing a coefficient of said polynomial, and solving said polynomial to thereby produce lens heights each corresponding to a particular image height input.
 40. The method as claimed in claim 39, wherein said corrected configuration calculating step comprises approximating the lens height distribution and a partial curvature distribution in a subscanning direction to a function model while storing a coefficient of said function model, solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce partial curvatures, and solving an arc function including said partial curvatures from positions in the subscanning direction to thereby produce heights.
 41. The method as claimed in claim 40, wherein said corrected configuration calculating step comprises: a modeling step of approximating the lens height distribution and the partial curvature distribution in the main scanning direction to a function model while storing a coefficient of said function model, and integrating said function model two times while storing a coefficient of an integrated function model; and a configuration calculating step of solving the function model at a plurality of positions of the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 42. The method as claimed in claim 38, wherein said sensitivity calculating step comprises superposing configuration errors only on the particular optical function surface by a simulation, determining, at a plurality of deflection angles, image heights, shifts of the focal point in the direction of an optical axis, and partial curvatures in a portion of said particular optical function surface through which a luminous flux passes, dividing said shifts of a focal point by said partial curvatures to thereby determine sensitivity of said shifts for a unit curvature, approximating an image height distribution and a sensitivity distribution to a polynomial or similar function model while storing a coefficient of said polynomial, and solving said polynomial to thereby produce sensitivity corresponding to an image height input.
 43. The method as claimed in claim 38, wherein said corrected configuration calculating step comprises approximating a lens height distribution and a partial curvature distribution in a subscanning direction to a polynomial: $C_{s} = {\sum\limits_{i = 0}^{n}{C_{si}x^{i}}}$

and using a coefficient C_(si) of said polynomial to solve an equation: ${z\left( {x,y} \right)} = {\sum\limits_{i = 0}^{n}{\frac{C_{si}}{2}x^{i}y^{2}}}$

for thereby calculating a corrected configuration at any desired point.
 44. The method as claimed in claim 38, wherein said corrected configuration calculating step comprises: a modeling step of executing numerical integration with a lens height distribution and a partial curvature distribution in a main scanning direction two times, and approximating results of sa d numerical integration to a function model while storing a coefficient of said function model; and a configuration calculating step of solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 45. The method as claimed in claim 37, wherein said lens height calculating step comprises determining by an optical simulation lens heights in a main scanning direction, at each of which a chief ray deflected by a particular angle is incident to a particular optical function surface, and image heights in said main scanning direction at each of which said chief ray is incident to an image surface, approximating a lens height distribution and an image height distribution to a polynomial or similar model while storing a coefficient of said polynomial, and solving said polynomial to thereby produce lens heights each corresponding to a particular image height input.
 46. The method as claimed in claim 45, wherein said sensitivity calculating step comprises superposing configuration errors only on the particular optical function surface by a simulation, determining, at a plurality of deflection angles, image heights, shifts of the focal point in the direction of an optical axis, and partial curvatures in a portion of said particular optical function surface through which a luminous flux passes, dividing said shifts of a focal point by said partial curvatures to thereby determine sensitivity of said shifts for a unit curvature, approximating an image height distribution and a sensitivity distribution to a polynomial or similar function model while storing a coefficient of said polynomial, and solving said polynomial to thereby produce sensitivity corresponding to an image height input.
 47. The method as claimed in claim 46, wherein said corrected configuration calculating step comprises approximating the lens height distribution and a partial curvature distribution in a subscanning direction to a function model while storing a coefficient of said function model, solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce partial curvatures, and solving an arc function including said partial curvatures from positions in the subscanning direction to thereby produce heights.
 48. The method as claimed in claim 37, wherein said corrected configuration calculating step comprises approximating a lens height distribution and a partial curvature distribution in a subscanning direction to a polynomial: $C_{s} = {\sum\limits_{i = 0}^{n}{C_{si}x^{i}}}$

and using a coefficient C_(si) of said polynomial to solve an equation: ${z\left( {x,y} \right)} = {\sum\limits_{i = 0}^{n}{\frac{C_{si}}{2}x^{i}y^{2}}}$

for thereby calculating a corrected configuration at any desired point.
 49. The method as claimed in claim 37, wherein said corrected configuration calculating step comprises: a modeling step of approximating the lens height distribution and the partial curvature distribution in the main scanning direction to a function model while storing a coefficient of said function model, and integrating said function model two times while storing a coefficient of an integrated function model; and a configuration calculating step of solving the function model at a plurality of positions of the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 50. The method as claimed in claim 37, wherein said corrected configuration calculating step comprises: a modeling step of executing numerical integration with a lens height distribution and a partial curvature distribution in a main scanning direction two times, and approximating results of said numerical integration to a function model while storing a coefficient of said function model; and a configuration calculating step of solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 51. The method as claimed in claim 37, wherein said sensitivity calculating step comprises superposing configuration errors only on the particular optical function surface by a simulation, determining, at a plurality of deflection angles, image heights, shifts of the focal point in the direction of an optical axis, and partial curvatures in a portion of said particular optical function surface through which a luminous flux passes, dividing said shifts of a focal point by said partial curvatures to thereby determine sensitivity of said shifts for a unit curvature, approximating an image height distribution and a sensitivity distribution to a polynomial or similar function model while storing a coefficient of said polynomial, and solving said polynomial to thereby produce sensitivity corresponding to an image height input.
 52. The method as claimed in claim 51, wherein said corrected configuration calculating step comprises approximating the lens height distribution and a partial curvature distribution in a subscanning direction to a function model while storing a coefficient of said function model, solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce partial curvatures, and solving an arc function including said partial curvatures from positions in the subscanning direction to thereby produce heights.
 53. The method as claimed in claim 51, wherein said corrected configuration calculating step comprises approximating a lens height distribution and a partial curvature distribution in a subscanning direction to a polynomial: $C_{s} = {\sum\limits_{i = 0}^{n}{C_{si}x^{i}}}$

and using a coefficient C_(si) of said polynomial to solve an equation: ${z\left( {x,y} \right)} = {\sum\limits_{i = 0}^{n}{\frac{C_{si}}{2}x^{i}y^{2}}}$

for thereby calculating a corrected configuration at any desired point.
 54. The method as claimed in claim 51, wherein said corrected configuration calculating step comprises: a modeling step of approximating the lens height distribution and the partial curvature distribution in the main scanning direction to a function model while storing a coefficient of said function model, and integrating said function model two times while storing a coefficient of an integrated function model; and a configuration calculating step of solving the function model at a plurality of positions of the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 55. The method as claimed in claim 51, wherein said corrected configuration calculating step comprises: a modeling step of executing numerical integration with a lens height distribution and a partial curvature distribution in a main scanning direction two times, and approximating results of said numerical integration to a function model while storing a coefficient of said function model; and a configuration calculating step of solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 56. The method as claimed in claim 37, wherein said corrected configuration calculating step comprises approximating the lens height distribution and a partial curvature distribution in a subscanning direction to a function model while storing a coefficient of said function model, solving said function model at a plurality of points on the particular optical function surface from positions in the main scanning direction to thereby produce partial curvatures, and solving an arc function including said partial curvatures from positions in the subscanning direction to thereby produce heights.
 57. The method as claimed in claim 56, wherein said corrected configuration calculating step comprises: a modeling step of approximating the lens height distribution and the partial curvature distribution in the main scanning direction to a function model while storing a coefficient of said function model, and integrating said function model two times while storing a coefficient of an integrated function model; and a configuration calculating step of solving the function model at a plurality of positions of the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 58. The method as claimed in claim 56, wherein said corrected configuration calculating step comprises: a modeling step of executing numerical integration with a lens height distribution and a partial curvature distribution in a main scanning direction two times, and approximating results of said numerical integration to a function model while storing a coefficient of said function model; and a configuration calculating step of solving said function model at a plurality of points on the particular optical function surface from positions
 59. The method as claimed in claim 37, wherein said corrected configuration calculating step comprises: a modeling step of approximating the lens height distribution and the partial curvature distribution in the main scanning direction to a function model while storing a coefficient of said function model, and integrating said function model two times while storing a coefficient of an integrated function model; and a configuration calculating step of solving the function model at a plurality of positions of the particular optical function surface from positions in the main scanning direction to thereby produce heights.
 60. The method as claimed in claim 37, wherein said corrected configuration calculating step comprises: a modeling step of executing numerical integration with a lens height distribution and a partial curvature distribution in a main scanning direction two times, and approximating results of said numerical integration to a function model while storing a coefficient of said function model; and a configuration calculating step of solving said function model at a plurality of points on the particular optical function surface from positions
 61. An optical element produced by a method comprising: a shift measuring step of mounting said optical element molded by infection molding to an optical scanning device, which is equivalent to an optical scanning device to be actually used, and measuring shifts of a focal point in a direction of optical axis at a plurality of image heights on an image surface; a lens height calculating step of calculating positions on a particular optical function surface that correspond one-to-one to the plurality of image heights; a sensitivity calculating step of calculating sensitivity in terms of proportional coefficients of partial curvatures and shifts of focal point on the particular optical function surface that correspond one-to-one to the plurality of image heights; a partial curvature calculating step of calculating partial curvatures on the particular optical function by using the shifts of the focal point and the sensitivity; a corrected configuration calculating step of outputting, based on the lens heights and the partial curvatures, a corrected configuration of the particular function surface in a form of coefficients of a function model or map data; a correction machining step of machining for correction a specular surface frame of a mold corresponding to the particular optical function surface; and a molding step of molding the optical element by using the specular surface frame corrected.
 62. The optical element as claimed in claim 61, wherein at least one optical function surface is substantially flat and constitutes the particular optical function surface.
 63. In an optical scanning device including an optical element, said optical element is produced by a method comprising: a shift measuring step of mounting said optical element molded by injection molding to an optical scanning device, which is equivalent to an optical scanning device to be actually used, and measuring shifts of a focal point in a direction of optical axis at a plurality of image heights on an image surface; a lens height calculating step of calculating positions on a particular optical function surface that correspond one-to-one to the plurality of image heights; a sensitivity calculating step of calculating sensitivity in terms of proportional coefficients of partial curvatures and shifts of focal point on the particular optical function surface that correspond one-to-one to the plurality of image heights; a partial curvature calculating step of calculating partial curvatures on the particular optical function by using the shifts of a focal point and the sensitivity; a corrected configuration calculating step of outputting, based on the lens heights and the partial curvatures, a corrected configuration of the particular function surface in a form of coefficients of a function model or map data; a correction machining step of machining for correction a specular surface frame of a mold corresponding to the particular optical function surface; and a molding step of molding the optical element by using the specular surface frame corrected.
 64. The optical element as claimed in claim 63, wherein at least one optical function surface is substantially flat and constitutes the particular optical function surface.
 65. A program for generating NC (Numerical Control) data to be used for machining with an NC machine a specular surface frame corresponding to an optical function surface of a mold for plastic injection molding, said program comprising: lens height calculating means for calculating, based on input data representative of shifts of a focal point measured at a plurality of image heights on an image surface of an optical scanning device, a lens height on a particular optical function surface corresponding to any one of said image heights; sensitivity calculating means for calculating sensitivity in terms of a proportional coefficient of a partial curvature and a shift of a focal point on the particular optical function plane that correspond to any one of the image heights; partial curvature calculating means for calculating a partial curvature by using the shift of a focal point and the sensitivity; and a corrected configuration calculating means for calculating, based on a lens height distribution and a partial curvature distribution, NC data for machining a particular specular surface frame to be corrected. 